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5_j . Abstract 

A model of evolution of bipartite quantum state entanglement is studied. It involves the recently introduced 
quantum block spin-flip dynamics on a lattice. We find that for initially separable states the considered evolution 
leads, in general, to entangled states. We also present a complete characterization of two-point correlation 
£f~) , functions for that type of dynamics to confirm enhancement of quantum correlation for the considered system. 

cn ■ 

> ', 1 Introduction 

m . 

£Q Entanglement^] first noted by Schrodinger and von Neumann as the characterictic trait of quantum mechanics 
I ■ Q , ||(| and by Einstein, Podolsky, Rosen || as a quantum superposition of two distinct states of physical systems 
| with the great impact on our understanding of the notion of reality in the atomic scale, is undoubtedly an essential 
t-H ■ feature of quantum mechanics. That concept can be relatively easy treated for the hamiltonian type evolution. 
Namely, writing the full Hamiltonian of a composite system, one should specify the interaction part explicitly. We 
. want to emphasize that it is exactly that part of hamiltonian, which is responsible for the evolution of entanglement 
O |' in the following sense: it is impossible that all factorizable states remain factorizable during the interaction unless 
, 1 * . the interaction part of the hamiltonian is trivial. 

The question of the evolution of entanglement for quantum stochastic semigroups is much harder task. To 
explain that point in detail, we recall that in recent papers, B. Zegarlinski and one of us, proposed a general 
scheme for a quantization of stochastic dynamics which describe interacting classical particles jl(J-|jl3|. In that 
scheme, guided by classical theory (cf. p|), a general recipe for quantum stochastic dynamics of jump type and 
diffusive type was given (see also [fj],|§],[5|],@])- In particular, to define the infinitesimal generators of Markov 
semigroups, the theory of generalized conditional expectations (in the Accardi-Cechnini sense) in the framework of 



5h ' quantum L p -spaces was used. In that way, a general scheme to produce, describe and analyze dynamical systems 
i with evolution originating from quantization of stochastic processes and such that their equilibrium states are 
prescribed (quantum) Gibbs states, was established. 

Having such a general plan of quantization of stochastic dynamics one may pose a natural question about its 
nontriviality. Under this notion we understand, first of all, that infinitesimal generator of such dynamics is not 
a function of the hamiltonian which can be extracted from the given Gibbs state. This requirement arises in 
a natural way from the methodology of construction of quantum stochastic dynamics sketched in the preceding 
paragraph. In fact, it has been shown that generators defined within the presented L p -space setting satisfy 
the above specification. On the other hand, the genuine quantum map should produce quantum correlations 
besides being well defined in non-commutative structures. The important point to note here is the idea of quantum 
correlations. Recently, the concept of coefficient of quantum correlations was introduced |]l6| and it was shown 
that such coefficient is not equal to zero only for non-separable states. In other words the entanglement is closely 
related to non-classical correlations between two subsystems of a composite system and such correlations arise from 
nontrivial interactions between the subsystems. Therefore, to have nontrivial quantum dynamics one should show 
that the dynamical maps under considerations are able to increase entanglement. 

The main difficulty in showing an increment of quantum correlations stems from the fact that the explicit 
form of interactions responsible for the transition rates used in the definition of quantum Markov generator is not 
known. Consequently, the present case is much more difficult from the hamiltonian one. To start a general analysis 
of evolution of entanglement for quantum stochastic dynamics, we begin with the block-spin flip type dynamics, 

x We refer the reader not familiar with the notion of entanglement to Appendix A, where some basic facts and definitions concerning 
separable and entangled states has been collected 



which is the main object of interest in this paper. Let us add that this dynamics can serve as a paradigm for a 
quantization of Glauber dynamics. Thus, we will be concerned with evolution of entanglement in concrete finite 
dimensional model of quantum block-spin flip dynamics. 

In order to carry out the analysis of entanglement we use two different approaches. Having fixed notations 
and given preliminaries (Section II), in Section III we consider an explicit example of low-dimensional jump- type 
system, showing that the block spin flip dynamics leads to the entanglement of the initial separable quantum state. 
In Section IV we study some properly chosen correlation functions and show that an enhancement of quantum 
correlation is typical for the considered dynamics (see Proposition 4.5). Finally, in Section V we will comment the 
obtained results. 



2 Quantum block spin-flip dynamics 

In the general approach to quantum jump-type dynamics of quantum systems on a lattice it is convenient to 
consider firstly the finite volume case - a system associated with a finite region A, and then to perform the 
thermodynamic limit with A going to Z d , where Z d denotes d-dimensional lattice (for mathematical details of 
algebraic description of quantum statistical mechanics see 0). However, as we have mentioned, the main scope 
of the paper is the analysis of evolution of entanglement and increment of entanglement is taken as a signature of 
nontrivial interactions between two subsystems. This phenomenon should be present both for finite and infinite 
subsystems. To simplify our exposition as much as possible we restrict to the essential ingredients of the finite 
volume case (for a general description see |l(J , [O] , jl3| ) . Thus, we shall consider a composite system I +11 associated 
with a region A = Aj U A]j, where A.;, i = I, II are disjoint finite subregions of the lattice Z d . To have a concrete 
dynamical system we will describe the construction of the block spin- flip dynamics related to the region A. To this 
end, we associate with A/ (A//) the finite dimensional Hilbert space 7Ya 7 = Hi {Hajj = H.2) as the space of its 
pure states, the set of density matrices S\ (S2) as the space of all mixed states, and the set of all bounded linear 
operators B(H.\) (8(0.2)) as the algebras of observables. Thus, the composite system A is described by Hi <g) H2, 
Sx ® S 2 and B(Hi) <g> B(H 2 ) = B(Hi <g> H 2 ), respectively. 

The reference state, playing the crucial role for classical and quantum case, is given here by the Gibbs state, 
i.e. with each region A (A/, A//, respectively) we associate the Hamiltonian H\ (H\j , H\ TI , respectively) and 

e -PH A 

subsequently the corresponding Gibbs state wa(-) = Trp\(-) with p\ = -ttt- = p, etc. 

Tre ^ 

Guided by the classical theory (cf. |S[]), where conditional expectations serve for the construction of jump type 
stochastic processes, we use their non-commutative generalizations (in the Accardi-Cechini sense) to define the 
infinitesimal generator £a,A/ = £ of the corresponding quantum spin-flip semigroup where the "block spin flip" is 
carried out on Aj (C A). Therefore let us introduce a map £a,A/ (= E) : B(Hi ® H2) ~^ B(H\ ® H2) defined as 
follows: 



E{A) = Tr x (7M7) (1) 

where 

1 = p^(Tr lP y h - 

with Tr\ denoting the partial trace (over the Hilbert space H\). Let us remark, that for infinite region case, 7 
is defined as analytic extension of non-commutative Radon-Nikodym cocycles (cf. [[l2|). Using the above defined 
generalized conditional expectation E we can introduce the following operator C defined on B(TL\ <g> H2) by 

C(A) = E(A) - A 

for A S B(Hi <8> 7^2)- The important point to note here is the form of the infinitesimal Markov generator C: there 
is no explicit term describing the interactions between regions A/ and An and this is the origin of difficulties in the 
presented analysis. Given a state w p , defined by a density matrix p, u> p — Tr(p(-)), one can define on B(H\ ® H.2) 
the following scalar product 

<<A5>>^Tr(^A*p*I?) 

Then, one can verify that (B(Hi ® H2), ((■, -))w p ) is the non-commutative Hilbert space (which can be called the 
quantum Liouville space). It will be denoted by L2(B(Hi ®H2),p)- Moreover, one can show that £ is a well 
defined bounded Markov generator such that 

((£(A),B)) u ^ = {{A,C(B))) u ^ 

It easily follows that the following semigroup T^(= T t ) — exp(t£) is a well defined Markov uniformly continuous 
semigroup such that it is self-adjoint on the non-commutative Hilbert space, the state lo p is invariant (with respect 



to T t ) and T t can be represented as the sum of the following convergent series: 

t 2 

1 + t£ + —£ 2 + •" 

with X being the identity operator. 

3 Evolution of entanglement 

In the next two sections we will look more closely at the time evolution of quantum correlations. Here, we 
will present the simplest nontrivial example clearly showing that quantum dynamics can produce this type of 
correlations. We will analyze a 2 x 2 system with block spin-flip dynamics. Thus, Hi = C — H2 and the Hilbert 
space of the composite system is given by C 2 ® (D 2 « C 4 . Let be an orthonormal basis in C 2 . We define: 

xx = ^=(&<8£i+6<8>&) 

%3 = 6 <8> a (2) 
a?4=^=(£i<8£i-6<8>$0 

One can easily check that {xi}f =1 forms the orthonormal basis in (D 4 . Let us define a faithful density matrix p 
on C 4 which is given by the formula: 

4 4 

P ■= Aj|a;i)(xi| A; > 0, J^ i = 1 ( 3 <* 
i=l i=l 

The just defined p will play the role of the reference state (cf. Introduction). We will need the following the well 
known fact which can be verified by straightforward calculation. 

Proposition 3.1 LetHi,H2 be Hilbert spaces. For every x,v € Hi andy,z£H2 we have: Tr\ (|a; ® y)(v <g> z\) — 
(x, v)\y)(z\ } where Tt\ is the partial trace over TL\ in Hi <8> TL2 and (•, •) denotes the scalar product (in Hi). 

In the sequel we shall identify Tri(-) with its embedding into Hi <S> H2, defined as 1 <E> TVi(-). Now, let us consider 
the quantum spin- flip type dynamics T t for our model. Its infinitesimal generator is defined as (cf. Section II): 

£(/) = E(f) - f 

where / G B(Hi <S> H2), while E was defined in the previous section. 

The dynamical semigroup T t has the additional property: as A C £2 (-4, p), we have T t A C A (Feller property), 
where A = B(Hi ® H2), ^{A, p) is the non-commutative Hilbert space. 

It follows that the Feller property allows to study the following duality problem: We may consider the time 
evolution T t as the family of maps T t : A — > A, then we can apply the standard equivalence between Schrodinger 
and Heisenberg picture to determine the evolution Tf of a state a. To this end we define: 

Tr(Tf(a)f) := Tr(aT t (f)) 

for any state a and any observable /. Therefore, we are able to describe explicitly the time evolution of states for 
that type of dynamics, which is given by the following mapping: 



a -> T?{o-)=o + t\E«{o)-o-J +■■■ (4) 

where the dual E d of the infinitesimal generator E is defined by the equality: 

Tr(E d (a)f)=Tr(*E(f)) 
and the series in the right hand side of (H|) is convergent. Using ([j]) we can write: 

Tr(E d ia)f) = Tr{aT n {Y fj)) = Tr(Tri(a) 7 */7) - Tr (7^1(^)7*/) 

Thus, we get: 

E d {a) = -<Tri{ah* = p? (Tnpy^no-tTnpy^ p? (5) 



Let us put a — a ®a, and recall that p was defined by M). Obviously 



Tr x o = lj® a 11 



Using Proposition 3J_ we can easily calculate Tr\{p) = J2t=i ^%Tri\Xi) {xi\. We have 



IYi|a;i) (aui | = -Tr t |6 ® 6 + 6 ® 6) <f l ® & + 6 ® 61 = 

i [Tri|6 ® 6><6 O £i| + Tn |6 ® 6><6 ® 6|+ 
+ Tn |6 <8 6) (6 ® I + Tn |6 ® 6) (6 ® 61] = 

i[iz® (16X61 +6X61)] = 5 [in] 

For simplicity, we will denote 1/ and In briefly by 1 when no confusion can arise. Analogously: 

Tn\x 2 ){x 2 \ = Tn|a ® 6X6 ® 61 = 1® 16X61 
r ri I13) (s 3 1 - Tn |6 ® 6) (6 ® 6 1 = i ® 16) (6 1 

Tn |x 4 ) (x 4 1 = ^Tn |6 ® 6 - 6 ® 6) (6 ® 6 - 6 ® 6 1 = \ [i ® i] 



Eventually, we obtain: 



2 

We introduce the following notation: 



Tr lP = yl + A 2 1 ® |6)(6I + A 3 1 ® |6>(6I + y 1 = 



A ' + A4 + \s) i ® 16X6 1 + ( + A 2 1 1 ® |6> (61 



2 

Inserting (|), @, and (@) into (|) we get: 



Xi = 7, 1" A 3 X2 = h A 2 



£<»= ^A?!^)^! • 1® [xfl6X6l+xII6X6l] 



• (iscr") • (^1® [x?|6X6l + xII6X6l]J • (^A?^ 

Now, suppose that 

°" J/ = «I6X6I + &I6X6I with a > 0,6 > 0,a + b= 1. Then, 

^V) = ( E A ? Nfcl ) ■ ( 1 ®«xil6X6l + i®&X2|6X6l) ■ ( 



x i=i 7 x i=i 

Using (0) and performing some lengthy calculation one can obtain: 



£ e V®( a l6X6l + &l6X6l)) 



aXi(AF+A|) 2 + 6 X2 (Af -A|) 2 



16X61® 16) (61+ 



+ (axi + b X 2) (Ai - A 4 ) |6) (6 1 ® 16) (6 1 + (0x1 + H2) (A x - A 4 ) |6) (6 1 ® 16) (6 1 + 



«Xi(A? - Af) 2 + 6 X 2 (A? + Af) 5 



16X61 ® 16X61 



+ axiA 3 |6)(6l ® 16X61 + &x 2 A 2 |6X6l ® 16) (6 1 

Now, we assume that A2 = A3 which implies Xi = X2 = X- Then, 



£> J ®(a|6X6l + &l6X6D) = x 



1 1 , 



[o(Af + A|)' + 6(Af - Af )'] |6X6I ® 16X61+ 



(Ai - A 4 )|6)(6I ® 16) (6 1 + (Ai - A 4 ) |6)(6l ® 16) (61+ 



1 1 , 



+ [a(A? - Af ^ + 6(A? + A|) z ] |6)(6I ® |6)(6I + 

aA 2 |6)(6l ® 16X61 + &A 2 |6X6I ® 16X61 

Performing some easy but tedious calculations we arrive at the following decomposition: 

4 

E d {a) = Y,Mvi)(Vi\ 
»=i 

with 

Ai = xA+ A 2 = xb^2 A 3 = x«A 2 A 4 = xA_ 



where 



and 



X± = A + C±X_ X = V{A _ C) 2 + AB 2 



A;=a{\\ +\\f+b(\\ -XI) 2 S:=Ai-A 4 C := a(A 1 i - Af) 2 + b(Af + A| ) 2 



while {yj}f = i is the orthonormal basis defined as below: 



yi = ??+6 ® 6 + «+6 ® 6 
y 2 = 6 ® 6 
2/3 = 6® 6 

2/4 = »7-6 ® 6 + «-6 ® 6 



where 



V2B 



\/X 2 T (A- C)X 



k± := 



-(A-C)±X 
W* 2 T (A- C)X 



The above decomposition of E d (a) is well defined for Ai > A 4 . In particular, we have A± > 0, rj±, «± =^ 0. 

Now, we are in position to examine the separability of the state E 4 ^ 1 ® ( fl |6)(6l + ^16) (6 1))- We will use the 

simple argument presented in H. Define: 



We observe 



and 



Let us put 



£ d (^)=Ai|yi)(2/i| + A4|j/4)(t/4| 
E$(a)y 2 = = E*(a)y 3 
(i$(cx)6®6,6®6) = Air?+K+ + A 4 r/_ K _ =A 



where A, (Bj) are positive operators in B(Hi) (£>(W 2 )). Then, 

\\E$(a) - S\\ > A J2W A U^\\ \\ B h*W \\ A J&\\ W B Ui\\ 

3 

Hence, 

\\E*(*)-S\\>±>0 

Consequently, E^(a) is not an element of the closure of separable states. To show that (1 — t)a + tE d (a) is an 
entangled state we start with recalling some well known facts from the theory of partially ordered topological spaces. 
Let us consider A sa = {» € A : a = a*}, where A — B{TL\®TL2) — B{TLi)®B{TL2) as a real finite dimensional Banach 
space {Hi and 7i 2 were assumed to be finite dimensional Hilbert spaces). The set {conv(B+(Hi) <£> B+(H.2))} closure 
will be denoted by V. Clearly, V is a proper generating cone in A. Further, by general argument (or by direct 
proof) one can note that intV ^ <f> (where int stands for interior of the set). Define V\ = V fl K(0, 1), where 
7^(0, 1) = {x G A : ||a;|| < 1}. One can verify that V\ is convex, compact subset with non-empty interior. 
Therefore, it is homeomorphic to unit ball, while its boundary dV\ is homeomorphic to unital sphere. Thus, it 
follows that, in general, a convex combination of p x e V\ and p 2 V\ is not in V\. By applying the above facts 
to the case p\ = (1 — t)a, p 2 = tE d (a), and taking into account the equivalence of the trace and operator norms 
in finite dimensional case, one can draw the conclusion that within the perturbation calculus in the first order a 
separable state a evolves to the entangled state (1 — t)cr + tE d (a). 



4 Factorization of two-points correlation functions 

In this section we will proceed with analysis of factorization of certain correlation functions. To clarify the relation 
of that topic to evolution of entanglement we start with establishing the relation between factorization and existence 
of non-quantum correlation between two subsystems A\ , Ai . Having this as well as the relation between quantum 
correlations and entanglement, we give a detailed analysis of factorization of two-point correlation functions. 

Let us begin with a classical case, i.e. let Ai, i = 1,2, be two abelian C*-algebras with identities. Consider 
Ad = A\ ® A 2 and let lj : A c i — > C be a state. The abelianess of Ai, i = 1, 2, implies 

Ad = C(fii) ® C(fi 2 ) = C(fii x 2 ), 

where fli, i — 1,2 are compact Hausdorff topological spaces. By Markov- Riesz theorem there exists a probability 
measure p on fli x fi 2 such that 

w(oi ® a 2 ) = p(ai2), 

with ai2 = i(a\ ® a 2 ) and i being an isomorphism between C(f2i) (£>C(£l2) and C(fl\ x Cl 2 )- Consider the truncated 
correlation (or equivalently the second Ursell function) 

Coj, ai ,a 2 = u(ai ® 1 • 1 ® a 2 ) - w(ai ® 1)^(1 ® a 2 ) 

for ai € -4i, i = 1,2. Due to the fact that each (classical) measure can be *-weakly approximated by finitely 
supported probability measures one has 

Cu, ai .a 2 = 1™ [Mn(ai ■ 02) - A*n(ai)Mn(tt2)] , 

n 

where for each n, p n = (n) , and c^m is a Dirac (point) measure on fii x f2 2 . Therefore S^(n) — 

5 ( n ,i) x <5 („, 2 ), where we have used the following notation: f^i x 2 3 = (wj™' 1 ^, w|™' 2 ^). Consequently 

Cu,a 1 ,a 2 = li m 
n 

Thus, each classical truncated correlation C Wi0li<l2 can be approximated by a difference of "separable states" and 
product of one-point functions. The important point to note here is that from the very beginning it is necessary 
to determine subalgebras (here Ai,A2) - then we can define correlations with respect to this fixed partition. We 
recall that we are studying block spin- flip dynamics, where the spin- flip is carried out over the region A/. Thus, 
the partition was fixed. Moreover, by definition, a separable state has the factorization property! 

Now, let Ai be arbitrary (non-commutative) C*-algebras with identities. Let w be a separable state on A = 
-4i<8>-4 2 , where in our case A = B{Hi)®B{H 2 ), Ai = B(Hi), i = 1,2. Then, guided by the above general observation 



,i=l 



concerning classical systems, we say that a separable state u> = \p\ ® Pi 1 encodes classical correlations with 
respect to the partition of A =< A\ <£> 1, 1 <8> A2 >, a C*-algebra generated by Ai (g> 1 and 1 ® .4 2 - 

Again, let w be a separable state with respect to the partition ,4=<.4i<g>l,l<g>.42 > and T t : A — > A be a 
spin-flip type dynamics. We wish to consider the truncated correlation (now for our quantum dynamical system) 

C q J T gJ =u(gT t (f))-u(g)w(T t (f)) 

with g £ id-Hi ® ^2, / G A\ ® id-H 2 . Throughout the rest of the paper we shall assume that observables / and g 
are of the form: 

/ = F <g> id H2 , c/ = <g> G 
with F G B(Hi) and G G 6(0.2), i.e. / (5) is an element of the subsystem I (II respectively). We observe: 

<W - w(ffTt(/)) - w(5V(r t (/)) - w(gf) + MgC(f)) + f ^(gC(C(f))) + ...- u>(g)u>(f) 
— ri -±ir L - x -i-* 2 n L > 2 < ^ r iL ' 2 j. 

where 

Cl aJ ^uj(gf)-uj(g)uj(f) 

(tii = <%£(/)) <%M£(/))> = Lo(gC(c(f))) u {g) u (C{C(f))), . . . 

We observe: 

i) G^ j measures the " classical" correlations between g and / with respect to the partition A\ , ^2, for time 
equal to 0. 

ii) We can not say the same for C^'g p C^'g f e * c -> smce w i s separable with respect to the partition Ai ® 1, 
1<S>.42, while g £ A2 but C(f) £ A\ We would have "classical" correlations if u was of the form uj = J2i ^i^f , 
where w{ is a state on a G*-subalgebra containing C(f) while ujf 1 is a state on A2, and similarly for higher order 
terms. In other words, the evolution (its infinitesimal generator) leads to a deviation from the original partition of 
the composite system and the new partition does not fit to our definition of the spin-flip operation over the region 
/. In particular, the given separable decomposition of u is not adapted to be a measure of classical correlations 
between g and £(/). 

Therefore, nonzero value of uj(G£(F)) = , can be taken as an indicator for the increment of non-classical 
correlations (with respect to the partition < A\ ® 1, 1 ® A2 >)■ 

To elaborate that point a little bit further let the state ui be given by the density matrix p of the form 

P = I>P?®pf. (8) 

i 

with p\ and p{ T being the states of the subsystem I and i7, respectively. Consider uo(GE(F)) — J2i KTr^pl ® 

p" ■ G ■ E(F)j with G £ A2, F £ Ai (with a natural embedding of Ai and A 2 into A). Assume G = HH* with 
[H > pj I }=0. Then, 

oj(GE(F)) = \Trp\ ® Hp\ I H*E(F) = ]T X^ 1 (HH*)Trpl ® pf ' H £(F) = 
= TrE d (^2 l \ i p I i I (HII*)p I i ® p\ I,H ^j -F 

i 

with pf = pi 1 / Pl (H*H). We know, by Section III, that E d ^Y Ji ^iP I i I (HH*)p I i <g> p[ 7,H ) is not, in general, a 
separable state. In other words, any reasonable factorization does not hold. In particular, 

E d ( ]T KpWHH^pl ® pf'") (F ® id„ a ) ^ £ Kp{ I (HH*) P i(F) 

i i 

We can summarize the above consideration as follows: 

An entangled state encodes non-classical correlations. Even in the simplest case of 2 x 2 system the spin-flip 
dynamics can lead to entangled state (cf. Section III). Moreover, in general, for the entangled state, the second 
Ursell function fails to have the factorization property. 

Before performing the promised analysis we want to make an additional observation, which justifies our as- 
sumption that it is enough to consider only observables from subalgcbras while considering correlation function 
lo(GE(F)). 



Remark 4.1 We note that the equality 



= w(GC(F)) = uj(GE(F)) - u(GF) (9) 

for all G G A and F G A2 is a very strong condition, since it implies E — id, so the spin-flip dynamics would be a 
trivial one. Therefore we restrict ourselves to study much weaker condition 

uj(GE(F)) ^ uj(GF) 

for G G Ai , F G Ai and a separable state ui on A\ ® A2 ■ 

Consequently, we shall analyze conditions under which the function < G,F >^ ui(GE(F)) does not factorize, i.e. 
@ does not hold for G G Ai, F G A2- Again, we assume that the state p is of the form (||). From now on we make 
the assumption that the density matrix p is an invertible one. This assumption stems from the general strategy of 
constructing quantum maps. Namely, we associate with the quantum system the quantum Hilbert space L%(A, p) 
with a given reference Gibbs state. Obviously, any Gibbs state has the assumed property. We also note that if p 
is an invertible density matrix, then Tr\(p) has this property too. This observation will be used throughout this 
section. 

As a result of longish (however not difficult) calculations we get the following characterization of the considered 
correlation functions: 

Let {tpi®4>j} be the orthonormal basis of Hi®H2, where {ifii} is arbitrary orthonormal basis of Hi while the basis 
{4>j} in H.2 is such that Tr\(p) is diagonal. Then we have: 

W)ff) P = EfE4^\/i)(^l F l^)(^|G|^) (10) 

where 



p 

klji x pq 



a j — Prrjj 
r 

and the matrix elements in the basis {tpi® are: 

Ppqrs := (<Pp ® fc\p\<fq ® 4>s) Ppqrs ■= (fp ® <l>r\p* \<Pq ® 0s) (H) 



Moreover, aj are the eigenvalues of Tr\ (p) . 



4.1 General characterization of factorization 

In order to examine the just described correlation functions, we assume as before that p, p 1 and p 11 (with or 
without indexes) denote states on Hi <8> H2, Hi and H2, respectively. 

It is evident from the definition of C^'j, g that the factorization of the correlation function (E(f)g} is not 
dependent on a particular choice of the decomposition of p. This allows us to examine the factorization for any 
decomposition of p with conclusions valid for any other decomposition. As the next step of mathematical framework 
for characterizing of correlation functions, we want to give the necessary and sufficient condition for the factorization 
of such functions. 

Proposition 4.1 Let p = ^ \p\ ® pj 1 an d f = F ® i&H 2 , g — idrti ® G, dimHi — n, dimH^ = m- Let {^i}™ = i 
be arbitrary orthonormal basis of Hi and {4>j}™=i be orthonormal basis ofH2 such that p = Trip is diagonal. Then 
the following conditions are equivalent: 

(i) correlation function can be factorized, i.e. 

(E(f)g) p =Y,H F ) p{ (G) pV 

i 

(ii) for every fc, / G {1, 2, . . . ,n} and j, i G {1, 2, . . . , m} the following equality holds: 



Plkji =2^1^ Ppk lq Plpqi\ " ( 12 ) 

p= l q=1 V 2^ V =lPvvii 

where pi kji := (ipi ® (j>j\p\tpk ® fc)- 

If the condition (ii) holds for some basis {fi} then it holds also for any other basis {&} (basis {<pj} unchanged!). 
Conversely, if (ii) does not hold for given basis {<Pi} then it does not hold for any other basis {^Pi}. 



Proof. Take the bases {<fik}, {4>j} such as described in the proposition. Calculate i (p) u m the basis 



{ip k ®(j>j}. Using (|ll|), we get: 

i i I j kljr 



Pvvjj 



From the general form of the considered correlation function (cf. (JlCj)) we have: 



kljr v pg 



Substituting _F = |y>fc)(y>i| and G = \4>j/\4 > i\ with k,l,j,i arbitrary, it is easy to verify that the right-hand sides of 
the last two expressions equal if and only if (|l2|) holds. 



4.2 Factorization and quasi-classicality 

It turns out that the sufficient conditions for factorization of correlation function can be connected to the 'quasi- 
classicality' of the considered state - the notion to be precised in the following definition. If p = J2 i Xipj ® pj 1 
is given decomposition of a separable density matrix, then p can be considered as classical if {pj} and {pj 1 } are 
abelian families of density matrices. Below, we define weaker conditions for families of density matrices, which are 
essential for the subsequent considerations. 

Definition 4.1 Let Ti be the Hilbert space, dimTL — n. Let {pi}fLi be a family of density matrices on Tt and 
{\}iLi be (strictly) positive numbers (Aj > 0), such that 53i=i ^» = !■ Define p := ^2n\pi- Let {</>j}" =1 - 
orthonormal basis consisting of eigenvectors of po and {cj}™ =1 - corresponding eigenvalues. We will say that the 
family {pi} is K- quasi- abelian (1 < K < n) if and only if the following condition holds: there exists the partition 
{Ap}p = i of {1, 2, . . . , n}, \A P \ > 1 (where \A\ denotes the cardinality of the set A), such that: 

(i) for any p £ {1, .... K} and r, s £ A p we have c r — c s and for any r £ A p , s G A q if p ^ q then c r ^ c s 

(ii) for any fc,i £ {!,...,«} and i f= {!,..., N} if Ck ^ ci then ((/)k\pi \<fil) — 



Remark 4.2 Definition ^A. can be rephrased in a simple way as follows. The family {pi} is K- quasi- abelian if and 
only if pi commutes with spectral projectors of po (so with po itself) for each i = 1, . . . ,N. However, we shall use 
properties (i) and (ii) in the sequel, so the original formulation is more convenient. 

It is an easy observation that if \A P \ — 1 for some p, then the corresponding vector <p p is the common eigenvector 
for all the matrices pi. This means that the number of sets A p in the partition {A p } with the cardinality one equals 
the number of common eigenvectors for all the matrices pi. In particular, if {p{\ is if-quasi-abelian for K = n, 
then it is abelian in the traditional sense. 

On the basis of the above definition we can formulate the sufficient conditions for factorization of the correlation 
function using the above type of commutativity properties of the families {pj} and {pf 1 }. 

Proposition 4.2 Let p — Xipj (dpi 1 be a separable density matrix and f — F(E>id-j-c 2 , g — id^^G, dirriHi = n, 
dimTi.2 = m. Then the following implication holds: 

There exists decomposition p = Y]j—i ^jPj ® Pj 1 such that one of the conditions is satisfied: 

(i) {Pj 1 } "is K- quasi- abelian (K < n) and {f>j} is abelian 

(ii) {Pj 1 } is abelian 



th 



en ( (E(f)g) p = J2i^i(F) p '(G) p n (factorization of the correlation function) j 



Proof. We will prove the implication assuming (i) . One can prove the statement under (ii) by similar reasoning. Let 
us first prove that ( |l2| ) holds. To this end let us take the bases {<fi} and {4>j} such that pj are diagonal in the basis 

{<fi} and Y) j ^jPj 1 = Tri(p) is diagonal in the basis {<pj}- Then, the matrix elements pikij — whenever I ^ k. 

i ii r — ~ 

The same is true for pf ki j. Note that also ^2 pq P p kj q Pi pq i\l a j I ' a i = for Z 7^ A: since for every p,i,j and q either 

Ppkjq = 01 flpqi = °- HenCe > fOT [¥> * © h ° lds - N ° W ' let 1 = k - Writing EpqplkjqPlpqi = (Eg PpkjqPlpqi) ' 

we see that it can differ from zero if and only if p = k(= I). Taking into account that due to our specific choice of 
the bases, the equality J2 q PkkjqPkkqi = Pkkji holds, we can write (|l|) for I — k in the form p kkji = pkkji\/Pjj/pii 
with pjj := (<f>j\p\(j)j}, p — J2 S Xspl 1 (= Tri(p)). Note that because of the choice of {<j>j}, elements pjj are the 



eigenvalues of p. By assumption, the family {pj } is K-quasi-abelian. This means that either pjj ~ pu ((|l2|) is 



then satisfied in an obvious way) or pjj ^ pa. In the latter case we have (<A,'|jo"|<fo) = (cf. Def. [O]), which 



implies pkkji = J2 S ^s{ l Pk\p I s \<Pk}{4 l j\p I s l \4'i) = and, of course, ( p"2| ) holds. We have shown that (i) implies 0(|i~2| 

By Proposition LI we have (E(f)g') — J] ■ A^i^—r^G)— 7 and thus, by equivalence of the decompositions of p, 

p Pj Pj 

we get (E(f)g) p = J2 t X *{ F ) p i{ G ) p u which ends the proof. 



The implication in Proposition 4.2 can be partially inverted. We have 



Proposition 4.3 Let p = A^pf (8pf J } / = F®id-H 2 > 9 = id-n x ®G, dimHi — n, dimYii — m. Assume that there 
exists decomposition p = J^j ^jPj ® Pj 1 such that {f>j} is ahelian. Then the following conditions are equivalent: 

(i) For some decomposition p = ^jPj ® Pj 1 we have: 

{f>j} is ahelian and {pj 1 } * s K-quasi-abelian (K < m) 

(ii) {E{f)g) p = Y,M F ) P l{G) p{I 



Proof, ((i) (ii)) This implication follows from Proposition 4.2. 

((ii) => (i)) Let {(fii} be the orthonormal basis of 7i\, in which pj are diagonal and {(j>j} be the orthonormal 
basis consisting of eigenvectors of the reduced density matrix p — Tr\(p). Then there exists decomposition 
p = y\ XjPj (g) Pj 1 such that pj — \ ( fij)( ( fij\ (the family {pj} is abelian). The matrix elements of p are pkiuw = 
J^i^iViki&uw with inki = (ipk\pl\tpi) = S ik 5u and £ iuw = ((p u \pj \<f> w ). Obviously, {E(f)g) p = J2i^i( F ) p i( G ) p n 



implies (S(/)g r ) = Y^,j ^j( F )^r\G)^ii- By Proposition 4.1 it means that 



PW =Y,H PpkjqPlwi \ # ( 13 ) 
P =i q =i V Pti 

i 

where p^ is eigenvalue of p corresponding to (f>j. Since pkiji = if I =/= k (which implies p k iji — for / =/= k), 

i i i i 

@ reduces to p kk ji = E g =i PkkjgPkkgt y/Pjj/P«, and JC 9 PkkjgPkkqi = Pkkji- Hence, if @ holds, then for every 
j, « G {1, ... , m} such that j / i we have: if ,0^- 7^ p,, then pkkji = for k = 1, . . . ,n. But pfc/yi = means 
J2i hVlkkilji = 0. Since r) ik i = S ik 5 a , we have = ^Z; Mvikk&ji = J2i hSlk&ji = ^k£kji €kji = as from the 
definition of decomposition V/ A; > 0. Thus, the set {pi 1 } has the following property: if pjj ^ pa then for every 



fee {1, . . . , n} we have (ft k )ji = £kji — 0. From Definition LI it follows that the set {pj } is K-quasi-abelian for 



some K < m. 

4.3 Factorization and nondegeneracy of density matrix 

The next result provides a criterion for the factorization of correlation functions under the nondegeneracy condition 
specified below. To show this equivalence we need the following result: 

Lemma 4.1 Let p — Xipj <£> pj 1 , dirnHi — n, dimli.2 — m and the matrix elements p k lji satisfy in some product 
basis {ipi <g> 4>j} the following condition p k iji = whenever k ^ I (pklji — whenever j ^ i). Then there exists 
decomposition p = Y]j Xjpj (8 Pj 1 such that {f>j} is abelian ({Pj 1 } * s abelian). 



Proof. Let us consider the matrix representation of p, i.e. p — [Pkiji\ k[ _ 1 n - j i—i m (cjj) ) • Suppose that 
Pkiji — whenever k ^ I. Let A s := J^Lj p S spp- Consider the following decomposition of p: 

n 

P = Y J ~^P I s®P I s ( 14 ) 

s=l 

where 

-/ I w I \~in I T •• 1 D»«ji]«=i ifA s >0 

[0 if A s = 

One can easily check that ( |l4| ) is a well defined decomposition of p (in particular if A s = then p SS ji — for 
j,i = 1, • ■ ■ , m). Of course {p^} is abelian. The proof of the second statement is similar. 

Now we are in position to give the promised result which shows a relation between factorization of correlation 
function and the spectral properties (nondegeneracy) of density matrix. 



Proposition 4.4 Let p = \ip\ ® pj 1 , f — F <g> irf^, <? = id,U\ ® G, dirriH\ = n, dirnrii — m. Assume that 
the reduced density matrix p= Kpj 1 (= Tri(p)) has nondegenerated eigenvalues. Then the following conditions 
are equivalent: 

(i) there exists decomposition p = ^jPj ® Pj 7 smc/i £/iat {pj^} * s abelian 

j 

(ii) correlation function can be factorized, i.e. (E(f)g') p — J^i ^i(^) p i(G) p n 



PV 



Proof, ((i) (ii)) This implication follows from Proposition 4.2. 

((ii) => (i)) Let {<fi}™ =1 be arbitrary orthonormal basis of 7i\ and be orthonor mal basis of TL2 con- 

sisting of eigenvectors of p. If the correlation function factorizes, then from Proposition 4.1 we have pikji — 
1 1 _ _ 

S p =i Sq=i PpkjqPipqi V Pjj I Pa w ith /9,j - the eigenvalue of p corresponding to </>,-. This equality and self-adjointness 
of p imply: 



E4 




EI I 
2 2 
PrlisPk 




r - : I / - ^ ] P Wis P&; 

Pjj 




2 j PrkjsPlr 




Multiplying the above equality by ( ^M- I gives: 



PpkjqPlpq 

pq 




T.Prk 3 sPl 1 



, Pjj _,. . / Pjj 

p« V p 



Plkji Plkji(Pjj — Pa) — 



— Plkji 



-Plkji 



By the assumption pjj ^= pa, hence pikji = whenever j ^= i. Our Proposition follows then from Lemma 4.1 



The results of this section provide a natural and intrinsic characterization of the two-points correlation function 
for block spin-flip dynamics and for the initial separable state. But one question still unanswered is whether the 
factorization or non- factorization of such functions is a genuine property for the considered dynamics. To answer 
this question we want to show that there are a lot of separable density matrices for which the correlation function 
{E(f)g) can not be factorized. Namely, we have the following: 

Proposition 4.5 The set of density matrices such that the equality (E(f)g} = \\Fj i\G} u does not hold, 

Trip is invertible}. 



is dense in S sep where S sep = {p € S sep 



The proof of Proposition L5 is given in Appendix B. 

We want to complete this section with the observation that there exists a strict connection between the problem 
of factorization of the correlation function and the separability of the square root p? . Namely: 



Proposition 4.6 Let p 

conditions is satisfied: 



then p? is separable. 



ji 



If there exists decomposition Y]. Ajp^i 

CO iPj} * s a belian 
(ii) {p 1 /} is abelian 



p . such that one of the following 



Proof. Suppose that there exists decomposition J^j ^jPj ® Pj 1 such that {pj} is abelian (the proof for the case 
when (ii) holds is similar). Let {if i} be the orthonormal basis of Hi, in which pjj are diagonal. Then, there exists 
decomposition p — ^2 k Xkpj ® pj 1 such that pj = \<Pk)(y>k\- Define matrix p as follows: 



p--=l^Kpl®{pi 



n\ 2 



k 



Note that p is a linear combination with positive coefficients and matrices pj, and pj. are positive operators. To 
complete the proof we must show that p is the square root of p. We have: 



p-p 



E A l^(^) s 



1 



K 2 pi ® (pi 



II\2 



kl 



K^-piPi^iplTiplT 2 



E W • p 1 ^ {pk)*{pV) h = E ^ 



Note that the above sufficient conditions for the separability of the square root of p are essentially weaker than 
those for factorization of the correlation function. 



5 Conclusions 



Our analysis yields new information about the nature of quantum spin- flip dynamics. It was shown that this type of 
evolution can lead to entangled states, so to the family of states encoding quantum correlations. Furthermore, the 
detailed analysis of factorization property of the Ursell functions clearly shows that this is an expected phenomenon. 
It would be desirable to have the full description of evolution of entanglement but we have not been able to do 
this. The main difficulty in carrying out such a description is that we do not know the general characterization of 
interactions causing the spin-flip operation. However, we were able to give a detailed analysis of factorization of 
two-point correlation functions. We should also emphasize that in our case there is no point in considering such 
notions as decoherence. Although the region A/ over which we perform a local operation (e.g. a block-spin flip) 
can be macroscopic, we do not deal with the collective observables. This means that studying the evolution of 
entanglement as well as the correlation functions is the proper tool in investigating genuine microscopic properties 
of the considered quantum dynamics. These two approaches, studied in Section III and IV, respectively, are not 
equivalent. It is well known that if at least one of the two subsystems is a classical system (i.e. the underlying 
algebra of operators is commutative), there is no entanglement, even if the second subsystem is a purely quantum 
one On the other hand, studying the correlation functions allows us to "detect" quantum prop erti es of the 
block-spin flip dynamics also when either of the two subsystems is a classical one (cf. Proposition |4.3| ) . Thus, 
in that aspect, the correlations based approach is a more "subtle" tool as far as we examine the problem of the 
considered dynamics being the genuine quantum map or not. However, even the detailed analysis of the second 
Ursell function (the analysis of the third, fourth, ... Ursell functions is a very difficult task) can not supply enough 
information about the time evolution of states of the system. From this point of view the approach given in Section 
III is much more fruitful. 

Turning back to the block-spin flip dynamics, we pointed out that its Markov generator should encode coupling 
between the region Aj and its environment A]j. Indeed, the results of Section III and IV say that the considered 
dynamics leads to enhancement of correlations. This means that the effect caused by the block-spin flip operation 
is strong, and it leads to coupling between two subsystems, therefore to nontrivial interactions. This enables us 
to interpret our result as another evidence that L p -approach to quantum dynamics is working well in the sense 
that it leads to a fruitful recipe for explicit construction of interesting genuine quantum counterparts of classical 
dynamical maps. 

Finally, we would like to remark that our results have been obtained for a low-dimensional model. Therefore, the 
expected and described properties of block-spin flip type dynamics follow exclusively from the noncommutativity 
of the underlying algebra of operators and have nothing to do with any transition from a finite to an infinite 
model via thermodynamic limit. Furthermore, we would like to emphasize that the presented theory has a fairly 
straightforward generalization to the infinite dimensional case as well as to other quantum jump processes. 

Acknowledgement: It is a pleasure to thank M. Marciniak for inspiring discussions. The work of (W.A.M) 
has been supported by KBN grant PB/0273/PO/99/16. 



Appendix A 



For the convenience of the reader we recall the definition of separable and entangled states in the general setting 
of Hilbert spaces. The density matrix p (state) on the Hilbert space Tii ®>Ti.2 is called separable if it can be written 
or approximated (in the norm) by the density matrices (states) of the form: 

P = ^2PiPi ®Ph ^(•) = ^J»i(w i 1 ®w i a )(-)J 

where p, > 0, Y^iPi = 1j Pi are density matrices on H a , a = 1,2, and (ojj ® ojf)(A g) B) = w}(A) ■ wf(-B) = 
{Trp\A)-{Trp 2 l B)=Tr{p]®pl-A®B}. 

Definition 6.1 Let p be the separable state on Tii ®1i,2- We say that every finite sum of the form p — \p\ ® p\ 
is a decomposition of the state p iff 

(i) Ei *i = 1 Vj A 4 > 

(») Kp\ ®p1 = p 

The state which is not separable is called non-separable. Denote by S the set of all states on Hi <X> H.2- 
Definition 6.2 Non-separable states are called entangled states. The set of entangled states is defined by 

^entangled = S\S se p 

where S sep stands for the set of separable states. 



For a discussion of physical aspects of that concept see jl8| , JlTj , . 



Appendix B 



In this appendix we give the proof of Proposition 4.5. Let us introduce the following notation: S n d C S sep , 
p G S nc i if and only if the eigenvalues of p = Tr\(p) are not degenerated, S n f C S sep , p € S n j if and only if 
(i?(/)<?) can not be factorized, S n df '■— S n d H S n f. 

Lemma 7.1 Let dimTLi, dimH.2 < oo. Then the set S n d is dense in S sep in uniform topology (equivalently, it is 
dense in any operator topology). 

Proof. Let p G S sep and p — J2iLi ^iPi ® pV some decomposition of p. Let {$jYj=\ be the orthonormal basis of 
H.2 such that p = Tri(p) is diagonal. Denote by e^ the eigenvalue of p corresponding to eigenvector 0j. Of course, 
we have e^ := (<fo p|^»). Without loss of generality we can assume that only one eigenvalue is degenerated. In 
particular, we can assume that e\ — e-x- Let e > 0. We will show that there exists p G S nd such that \\p — p\\ < e. 
Take rj such that < r\ < ^min{e, \es — e%\, . . . , \e m — e i | } ■ Define: 



V,= 



i=l,...,N 



\ ■= Ai(l -rj), p\ ■= pI, P 



ji 



V, Pn+i 



1 



dimTL] 



idHi, Pn+i : = 



and 



N+l 

p := E \p\ 

i=l 



Note that p is a well defined density matrix on Hi <8> H.2- Moreover, the reduced density matrix Tri(p) = 
Silt 1 ^iPi 1 = P(l- 7 l)+ r l\<l ) i)(4>i \ has only nondegenerated eigenvalues e\ = ei, e 2 = e 2 (l-?7), . . . , e m = e m {l-rf), 
so p G S'nd. The lack of degeneracy stems from the choice of n because for all i = 3, . . . , m we have | — ei | > 2n 
and, evidently, \e.iT]\ < ?y. Now, suppose that ei = e^ for some j G {3, . . . , m}. We have: 

ei = e) ei = ej(l — rj) e\ — ej = —ejij => 



ei 



| e j- 77 1 2i] <r] <^ 7? < 



which yields a contradiction, since rj was assumed to be positive. 

To complete the proof we must check that the inequality \\p — p|| < e holds. Indeed: 



N 



N+l 

i 1 - E ^ 

i=i 

JV 



= E^^'-Et 1 - 7 ^ 



i=l 
JV 



JI 



VPn+i®Pn+i 



'Pi 1 ~V Pn+i® Pn+i 



= \\v^2^ipi 

i=l 

= v\\p~Pn+i ®Pn+i\ \ < v(\\p\ \ + \\Pn+i®Pn+i\\) < 2 V < e 



Lemma 7.2 Let dirriHi, dintH.2 < 00. TTien the set S n df is dense in S n d in uniform topology (equivalently, it is 
dense in any operator topology). 

Proof. Let p G S n d and dimH.2 = m. Suppose that (E(f)g}j can be factorized. Then, from the relation between 

factorization and nondegeneracy of density matrix (see section [Tj^) there exists decomposition p — J2iLi ^iPi ® pV 
such that {pf 7 } is abelian. We can assume that N equals dimTL2 and pf 1 — |<fo)(0i|, where {(f) j} is the orthonormal 
basis of Ti.2 such that p = Tr\(p) is abelian. Denote by e$ the eigenvalue of p corresponding to the eigenvector <fii 
(we have ej = (<^i|p </>«))• Without loss of generality we can assume that eigenvalues of p are ordered decreasingly, 
i.e. ei > e-x > . . . . 



Let e > 0. We will show that there exists p € S n df such that \\p — p]| < e. Take r\ such that < rj < e/2. Let 
{Vi}™=i ^ e arbitrary but fixed orthonormal basis of Hi- Define: 



i,...,n A, :— Aj(l 


- v), 




= pi, pi 1 


= pI' 


Pn+i-= 


W 1 ' 


(</>i| 


+ -\<P2)(V>2 





— 1 1 1 j 

Ajv+2 := 2^ PN+2 '■= \<P2)(<P2\, Pn+i-= + ^l^)^]- j\^>l)((P2\ 



<p2)(<Pl\ 



and 



JV+2 

i=l 



II 



Note that p is well defined density matrix on 7ii®H.2- Moreover, IYi(p) = Yl!i=i ^iPi 1 = ?(l — v) + 5 ? 7(| ( / , i)( < / , i | + 
\4>2){(t>2 1) has only nondegenerated eigenvalues ei = ei (l — -^77) , — e2(l — 517) , eb = 63(1 — 77), . . . , e m = e m (l — rj), 
sopeS nd . 

Now we aim at showing that p £ S n df- It is enough to show that there is no decomposition p = X)j=i Ajpf <8 p 7/ 
for which {pf 7 } is abelian, since due to the relation between factorization and nondegeneracy of density matrix 
(cf. section [4.3| ) it is equivalent to the fact that (-E(/)<?)~ does not factorize. Suppose that there exists such a 

decomposition with {pi 1 } abelian. We can assume that N equals dimHi and pj 1 — \4>i)(<j)i\. Then, we have: 



JY 



N+2 



2ll 



8=1 t=l 



i=l 



X (Xipj - Xipj) <g> j^X^I = Ajv+ijo^+i (g> p^f +1 + Ajv+2Pw+2 ® pfv^ 



Since (Ajpf — A^pf ) ® |^>i)(</>i| are linearly independent, we have \ip{ = \ip\ for i = 3 
following equality: 



which leads to the 



(\Pi - Ai/3f) ® |&}(<fc| = Xn+iPn +1 ® pfv+i + Aw+2Pw+2 ® PiV+2 



i=l 



Denote the left-hand side and the right-hand side of the above equality by L i P, respectively, we have: 

(<pi <g> 4>\\L\ip\ ® 0a) = and (<p! ® ^Ipltpi ® <^ 2 ) = ~»7 ^ 

8 

which yields a contradiction. Thus, p € S n( 2/. 

To complete the proof we must check that the inequality ||p — p\ \ < e holds. Indeed: 



N 



N+2 



x *p* ® pi 1 - x *pI ® # 7 



i=l 

AT 



i=l 
JV 



hp{ ® P- 7 ~ X (1 ~ J7)A,pf ® pf - ^77 plv+i 



i=l 
N 



1 PW+I ~ 2^ ® pfv+2 



= I A. ( pf ® p 77 - ir/pfv+i <g> pfv 7 +1 - -VPn+2 ® Pw+2 

i=l 

= »/| |P - ^PW+I ® P/V+l - ^PiV+2 ® PW+2I I 

< ^(l|p|| + d|p^ + i®P^ + ill + dl^+a®i3^+all) < 2 ^< e 



Proof (of Proposition jj-^l )- The following inclusions hold: S n df C S'nd C S'sep- According to Lemma 7.1 and Lemma 



7.2, S'nd is dense in S sep and S n df is dense in 5 n rf, respectively. It means that Sndf is dense in S sep . Moreover, we 
have: S n df C S n f C S'sep which implies that 5 ra / is dense in S'sep. The proof is completed. 
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